FUSION HIERARCHY AND FINITE-SIZE 
CORRECTIONS OF U q [sl{2)] INVARIANT 
VERTEX MODELS WITH OPEN BOUNDARIES 



The fused six-vertex models with open boundary conditions are studied. The 
Bethe ansatz solution given by Sklyanin has been generalized to the transfer ma- 
trices of the fused models. We have shown that the eigenvalues of transfer matrices 
satisfy a group of functional relations, which are the su(2) fusion rule held by the 
transfer matrices of the fused models. The fused transfer matrices form a com- 
muting family and also commute with the quantum group U q [sl{2)]. In the case of 
the parameter q h = — 1 (h = 4, 5, • • •) the functional relations in the limit of spec- 
tral parameter u — ► ioo are truncated. This shows that the su(2) fusion rule with 
finite level appears for the six vertex model with the open boundary conditions. 
We have solved the functional relations to obtain the finite-size corrections of the 
fused transfer matrices for low-lying excitations. From the corrections the central 
charges and conformal weights of underlying conformal field theory are extracted. 
To see different boundary conditions we also have studied the six-vertex model with 
a twisted boundary condition. 
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1 Introduction 



In statistical mechanics the commuting transfer matrix of two-dimensional lattice models 
is often used. This is the case because it shows obviously that the corresponding systems 
are integrable and can be solved exactly. The commutativity of the transfer matrices for 
periodic systems is easily derived from the Yang-Baxter equation Jl|, ||. Recently, the 
so-called open boundary condition has been introduced in the study of two-dimensional 
lattice models. To have the commuting family of transfer matrices for such systems we 
have to use the Yang-Baxter equation for the bulk and the reflection equation [[| for the 
boundaries. 

The exactly solvable models with non-periodic boundaries have been early studied in 
H) H 0) H Hi- After Sklyanin's work on the algebraic Bethe ansatz of six- vertex model 
with the open boundary condition j3|, there has been increasing interest in exploring 
two-dimensional lattice models or integrable quantum chains with the open boundary 
conditions jl2l-||2^|. Recently, the boundary cross unitary has been derived in []23f and a 



bootstrap approach has been developed to two dimensional integrable field theory with 
boundary in p3), |4| (see also |Tl|, |25], [26|] for related works). In a vertex operator 
approach has been used to solve the semi-infinite XX Z spin with a boundary magnetic 
field. 

The spin- 1/2 XX Z chain is the Hamiltonian limit of the transfer matrix of six- 
vertex model. Similarly, the higher spin XX Z chains are the Hamiltonian limits of the 



fused transfer matrices of six- vertex model |27| |28], [29|, [30| . This corresponding relation 
also works for the model with the open boundary conditions f|], j31| . As the fusion of 
the Boltzmann weights of six-vertex model |27j so the fusion procedure of the boundary 



matrices has been expressed in The fused transfer matrices with the open boundaries 
are well defined and form the commuting families of the model. They satisfy a group 
of functional relations, which can be shown by fusion procedure. For the Andrews, 



Baxter and Forrester critical solid-on-solid models |T(J it has shown that the functional 
relations are useful to find the Bethe ansatz solutions and the finite-size corrections of 
the fused transfer matrices |42[. But, for two dimensional vertex models it seems 



not the case. The functional relations of the fused transfer matrices of six-vertex model 



with periodic boundary are given in [38[ and have been shown that they are not closed 
f55fl . However, the situation is very different for the case of the six- vertex model with 
the open boundary conditions. In this paper we study the six-vertex model with the 
open boundary conditions. The fused transfer matrices of the model commute with the 
quantum group U g [sl(2)]. The eigenstates of the transfer matrices can then be classified 
according to the spin S z = j. We focus on the interesting case of q h = — 1 {h = 4, 5, • • •). 
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We show that the functional relations of the fused transfer matrices are truncated in braid 
limit. This important fact allows the functional relations to be solved in thermodynamic 
limit. We have found the finite-size corrections of the fused transfer matrices for the 
low-lying excitations, which are the spin S sectors above the ground state. Therefore 
the central charges and the conformal weights of underlying conformal field theory have 
been extracted from the finite-size corrections. The central charges are given by 

3p 6p 



p + 2 h(h-p) 
and the conformal weights are given by 



/i = 4,5,-.. (1.1) 



2 

h — (h — p)s ) — p 2 



u(p — v) , , 

+ — ^ £ s = l,3,---</i-l, (1.2) 



Ahp(h - p) 2p(p + 2 

where p — 1, 2, • • • is the fusion level and v is the integer determined by 

s- 1 



v 



V 



p. (1.3) 



Here the brackets \_x\ denotes the greatest integer less than or equal to x. For the case 
of p = 1 these conformal spectra coincide with the results of the spin-1/2 XX Z chain 
with the open boundary condition [|11[| . It is interesting to notice that for generic p these 
are the conformal spectra of the spin-p/2 XX Z chain with the open boundary condition. 
In a similar way we also have obtained the finite-size corrections of the fused transfer 
matrices of the six-vertex model with a twisted boundary condition. 

The outline of the paper is as follows. In subsection 1.1 and subsection 1.2 we explain 
the six-vertex model and the open boundary condition and recall the Bethe ansatz solu- 
tion of the U q [sl{2)] invariant transfer matrix of the six- vertex model. In section 2 the 
functional equations of the fused transfer matrices of the model are presented and the 
corresponding Bethe ansatz of the transfer matrices are constructed. In section 3 we carry 
out the procedure of calculating the finite-size corrections of the fused transfer matrices 
and extract the conformal spectra of underlying conformal field theories. In section 4 a 
brief discussion is given. Particularly, we discuss the non-?7 g [s/(2)] invariant six- vertex 
model with boundaries. In Appendices we describe the fused Boltzmann weights of the 
six-vertex model and show the functional equations directly by fusion. 



1.1 Six- vertex model 



The six-vertex model is one of solvable lattice models in two dimensional statistical 
mechanics to prove tractable |SB|, |34|, |2|. It is ice-type model and has six nonzero 
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Boltzmann weights. These Boltzmann weights form a four by four i?-matrix 

(a(z) \ 
b(z) c(z) 



R{z) 



where 






V o 



a{z) 
b{z) 
c{z) 





c{z) b(z) 
a(z) J 



zq — z 1 q 1 



z — z 



f M) 



1.5) 



q-q 



depending on a parameter q and a spectral parameter z G C. This i? matrix acting on 
C 2 <8> C 2 solves the Yang-Baxter relation 



R V2 {x/y)R lz {x/z)R 2 ' i {y/z 



23/ 



R 23 (y/z)R r6 (x/z)R 12 (x/y) . 



>13 



;i.6) 



For the twisted boundary condition the row-to-row transfer matrix acting on a " quan- 
tum" space C 2N = C 2 ® • ■ ■ C 2 is defined by 



trl g 



;i-7) 



where m is an integer. It becomes the transfer matrix with the periodic boundary con- 
dition for m = and for m > it is twisted. The trace is taken in the "classical" space 
C 2 and the monodromy matrix 



U(z) 



R c > N (zq-\) ■ ■ ■ R c ' 2 (zq-\) R c ^(zq-\ 



is a two by two matrix in the classical space C 2 denoted by c and with their elements 
acting on the quantum space C 2N . The Yang-Baxter equation (1.6) follows that the 
monodromy matrix satisfies the quadratic relation 



R 12 (x/y) U (x/z) U (y/z) = U (y/ z) U (x/ z)R 12 (x/y) 
and the transfer matrix forms a commuting family: 

[ T m {z) , T m (y) } = . 



;i.9) 



1.10) 



This means that T m (z) is a generating function of commuting operators in quantum 
mechanics of one dimensional chains. Specially, the Hamiltonian 



Hry 



n=l V Z 



of the quantum XX Z spin-| chain is contained within this family. These a*, a y n and o z n 
are Pauli matrices and a^r = of ± ia%. 



n n n " 



The six-vertex model with the open boundary is given by recalling the reflection 
equation ||, 

R 12 (x/y)Kl(x)R 12 (xy)K 2 ± (y) = K 2 ± (y)R 12 (xy)Kl(x)R 12 (x/y) . (1.12) 
Sklyanin has shown in || that the transfer matrix 

T(z) = tr K + {zq 1 l 2 )U{z)K_(zq- 1 l 2 )U-\z- 1 ) (1.13) 
forms a commuting family 

[T(z), T(y)]=0. (1.14) 
Therefore it presents an integrable system with the boundary described by the reflection 



matrices K±(z). Particularly, the matrices K±(z) take the form |L2 



K±{z)=[ Z l ji), (1.15) 



then the integrable system possesses the quantum algebra U q [sl(2)] symmetry, which 
means 

T{y)S ± -S ± T{y) = and T{y) S z - S z T{y) = (1.16) 
for any y. Here the operators 5' ± and S z of the quantum algebra U q [sl(2)] 

S* = \{a z l+ a z 2 + --- + a z N ) (1.17) 

N 

S± = ^gK+-+<- 1 )/2 (a ± /2 ) g -K +1 +-+^)/2 (i.is) 

n=l 

satisfy 

S Z S ± -S ± S Z = ±S ± and S + S~ - 5^5+ = [2S z ] q (1.19) 

where 

qX _ X 

[ x \q = zr • 

q-q 1 

The Hamiltonian corresponding to the integrable system is the open XX Z quantum 
spin-i chain [|] 

H = EV<< + i + «+i + i ^<< + i)+^T~(^ + ^) ( L2 °) 



n=l 
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The U q [sl(2)] invariance Q1.16Q implies the quantum algebra U q [sl(2)] is the "symmetric 
group" of the XX Z chain p5f 



H S ± - S ± H = and H S z — S z H = (1.21) 
For generic values of q the representations of U q [sl(2)] are known to be equivalent to 



the ordinary su(2) representations |51, [5^]. The representation theory becomes more 
complicated for the special case of q h — ±1 (see E2l K3l BR Bl, [55|, |5H] for details). 



1.2 Bethe ansatz solution 

The Bethe ansatz equations and the eigenvalues of the six- vertex model with the twisted 
or open boundary condition have been given using the algebraic Bethe ansatz 0, [| [| [H| . 
Set z = e m and q = e lA . We recall that the eigenvalues T m (z) or T(z) of the transfer 
matrices T m (z) or T(z). For the twisted boundary case the eigenvalues are given by 

TM = e -"sm~( M + lA)n Sin( "-'' t -, A) 

^ sui(m - %) 

+ sin> - §A) Sin( T^ A) (1-22) 
^ sm(w - v m ) 

and these t>i, t>2, • • • , % are the solutions of the Bethe ansatz equations 

sm N (v k + |A) _ 2imA it sin(t; fc - v t + A) 
sm JV (y k - I A) ^ sm(t) fe - — A) 

For the open boundary case the eigenvalues are given by 

t(u) = sin(2 ^ + A) sm^ju + ia) n sin(M 7 m ~ A l sin S M + A) 

sm(2u) ^ siii(m - t> m ) sm(w + u m ) 

+ Sin(2 "~ A) sin 2 ^ - IA) fi Sin(M 7 m + A | Sin j M + ^ + A) (1.24) 
sin(2w) 2 ^ sin(w - v m ) sin(w + u m ) 

and these v\, V2, ■ • • , vm satisfy the Bethe ansatz equations 

sin(?; m + |A) \ 2jV _ H sin(t; m - + A) sin(^ m + t> k + A) 
sin(u TO -|A)y ^ s\n{v m - v k - A) sin(v m + v k - A) ' 



2 Fused models and their functional equations 



The fusion models can be built up by fusion [27| from the six vertex model. Suppose 
that R(p tq )(u) represents the i?-matrix of fused vertex (see Appendix. A) and then the 
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relevant monodromy matrix is defined by 

U(p, q ){u) = R c (l p) {u)R c ^ p) {u)---,R c ^ N p) {u) , (2.1) 

where p and q are respectively the fusion levels for vertical direction and horizontal 
direction of the square lattice and p, q = 1, 2, • • •. With the periodic boundary condition 
the fused transfer matrices 

T<*«\u) = trU M (u) (2.2) 

commute 

' T M (u) , T {p ' q '\v)] =0 (2.3) 
for each fusion level p fixed and any z and y. These q and q' can stay in different levels. 



The transfer matrices satisfy the following functional relations |3q] 

T^\ u ) T fel) (u + q\) = T (M+1) (n) + /£_! T^ q - l \u) , (2.4) 

where T^' -* [u) = J, the identity matrix, and the ^-dependent function f p is generated 
from the antisymmetric fusion of the Boltzmann weights. These relations are the su(2) 
fusion rule. They mean the relationship among the eigenvalues of the fused transfer 
matrices. In other words, all eigenvalues T^'^ of the fused transfer matrices are deter- 
mined by the relations with the initial solution T^' 1 ' . In the following subsection the 
similar idea is used to the transfer matrices with open or twisted boundary condition. 
Let T^ p ' g \u) be the eigenvalues of the fused transfer matrices with the open boundary 
condition or twisted boundary condition. We prove the following theorems. 

Theorem 2.1 (su(2) Fusion Hierarchy) Let us define 

t (p,o) = l T (g) = T (P,<i)( u + k x) fP = fP( U + qX) 

f p (u) = + A)0(u+ \p\ + X)u 2 (u)(j)(u- \p\) (2.5) 

where uj\{u), ^(it) and 4>(u) are given by $U%) - ft2J$) and ( \2U3j ). Then the su(2 ) fusion 
hierarchy follows 

T M T (i) = T ( q +i) + fP_ lT i*-Q (2.6) 

forq = 1,2, 

Theorem 2.2 (su(2) TBA) If we define 

t°o = (2.7) 
tl = T^ +l) Tt l) / q j\f kl (2.8) 

then it follows that su(2) TBA equations 

^(l+tna+tr 1 ). (2.9) 
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2.1 Bethe ansatz for fused models 



The Bethe ansatz solutions ( |1.22| ) and ( |1.24| ) can be written in the form of 



T(u)Q(u) = uj x {u)(j){u + \\)Q(u - A) + u 2 {u)4>{u - \\)Q{u + A) (2.10) 

using Baxter's auxiliary matrix Q which commutes with the transfer matrix T(u). The 
eigenvalue Q of the auxiliary matrix Q is given by 



Q(u) 



M 

J\ {sin(w — v rn ) for twisted boundary 

m=l 

M 

J\ {sin(u — v m ) sin(w + v m )} for open boundary 

V m=l 



(2.11) 



and 4>{u) is given by 

<p(u) = sin Ar (n) (2.12) 

with 

N for twisted boundary 

Af={ y (2.13) 
I 2A^ for open boundary. 

The functions Q(u) and <f)(u) are not directly related to the boundary. The boundary 
terms come in the expression of eigenvalues through the factors uji(u) and u>2(u) 

uj\(u) = ^(w) 1 = e _lA for twisted boundary (2-14) 

sin(2w + A) 

uj\{u) = uj2\—u) = — r — for open boundary. (2.15) 

sm(2?i) 

Here the simple case m = 1 has been taken for the twisted boundary. 



The fusion procedure |27|, |37|] (see |3(| (IT], for related works) shows that the 



fused transfer matrices T^'^iu) can be constructed directly from the unfused ones by 
applying the fusion projectors to the "quantum" space (see Appendix. B). The procedure 
implies obviously that the eigenvalue T^ p,lS) (u) for the transfer matrix T^' 1 ^) has the 
following form 

T^iujQiu) = ui(u)(p(u + \p\)Q{u - A) + uj 2 {u)(t){u - \p\)Q{u + A) . (2.16) 

The boundary conditions take the position in the "classical" space and thus oj\{u) and 
u>2{u) are not effected by the fusion procedure in "quantum" space. The function Q{u) is 
however dependent on the fusion level p. For example, the ground state corresponds to 
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take M = pN/2. The Bethe ansatz equations determining all these v i, t>2, • ■ • are given 
by setting 

T [pA \v k ) = . (2.17) 



To show the su(2) fusion hierarchy let us use semi- standard Young tableaux 38, 39, 

Q(u + k\ + X) 



40], Define 

k 

rri = U2(u + k\)<p(u + k\ — \p\) 
H 

for a single Young tableau so that 

Ti Vl = H + H = E □ 



Q(w + fcA) 

" / . . \ , / , , , , ,Q(u + kX - X) 

u x {u + kX)(j){u + kX + \px y K ' 

Q[u + kX) 



(2.18) 



(2.19) 



For a general one-row Young tableau, the numbers must not decrease moving to the right 
along the row, e.g. 

n 

(2.20) 



1 


1 


2 


2 


2 



Such a Young tableau denotes the product of the five labeled boxes defined by ( [2.1 8|) 
where it is understood that the relative shifts in the spectral parameters are given by 



u + A\ 


u + 3A 


u + 2A 


it + A 


u 



(2.21) 

and the zero superscript gives the shift in the most right box. Filling the numbers 1 
and 2 in this five-box Young tableau according to the rule that the numbers must not 
decrease moving to the right along the row, we get six numbered Young tableaux. Then 
taking sum of these six Young tableaux with the correct spectral parameter shifts ( |2.21| ), 
we obtain the eigenvalues T^{u). 

By a similar way the eigenvalues of the fused row transfer matrix at level q can be 
written as 



T, 



(?) 



T^\u) 



E 



(2.22) 



where the number of terms in the sum is given by the dimension of the irreducible 
representations of su{2) 



dim(g) = (q + 1) 
9 



(2.23) 



For example, the fusion level q = p = 2 case gives the eigenvalues of transfer matrix of 
the 19-vertex model 



(2) 







H 


1 


+ 


1 


2 





+ 


2 


2 



(2.24) 



where it is understood that the relative shifts in the arguments are given by 

o 



(2.25) 



It is straightforward to show that set 

_Q 

= uji(u + X)(f)(u + \p\ + \)uj 2 {u)(j){u — \p\) 



F{u) 



then T (9) given by (|2.22|) satisfy ( |2.6|) . This leads the theorem [2.1| . The proof also shows 
that the fusion hierarchy is compatible with the su(2) fusion rule 



□ 



(2.26) 



5-1 



9 + 1 



To show theorem [2.2| let us consider the triple 

Inserting the fusion hierarchy into the terms in parentheses this equation gives new 
functional equations 

T^T[ q) = J] fk + T^T[ q - l) , (2.27) 
which corresponds to the following su(2) fusion rule 



<8> 



(2.28) 



9-1 



9 + 1 



Then it is easy to see the theorem |2.2| by rewriting the fusion rule according to the 
definition of t q (u). 

The functional equations (|2.6|) and fl2.9Q in form are the same as the su(2) functional 
equations of A-D-E models BP], HSJ and the dilute A-D-E models p7 . 
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2.2 Zeros and poles of eigenvalues 



The functional relations have been shown to be very useful to calculate the finite size 



corrections of the fused transfer matrices ||!| . To solve the fusion hierarchy (|2.6| ) and 



we need to know the distribution of zeros and poles of these transfer matrices 
and . For q = p these possess the physical strip of the model. Inside the strip the 
ground state eigenvalues T^ q > do not possess any zero apart from those which are imposed 
by the fusion of the Boltzmann weights and the boundary. The zeros contributed by the 
Boltzmann weights are of order Af and those by the boundary are only of order 1. They 
list them as follows. 

zero[T {p ' g) (w)] = for q < p (2.29) 

zerofT^O)] = [f-J- 1 {-kX - \ P X} N for q > p. (2.30) 

for the twisted boundary condition and 

zero[T M {u)] = \j£*{-k\-±\} for q < p (2.31) 

zevo[T^\u)} = ^l{-k\-\\}[fiJ-\-k\~\p\Y N iovq>p. (2.32) 

for the open boundary condition. The zeros and poles of t^ q ' are determined by (|2.8|) . So 
we have 

(I) q<p-l : 

zero[t (p ' 9) (n)] = 

q-l q-l 

pole[t ip ' q \u)] = |J {-k\ - \pX - \} N |J {\p\ - k\} N (2.33) 

k=0 k=0 



(II) q = p : 

zero[t (p ' p) («)] = {-ipA} 



N 



p—1 p— 1 

pole[t (p ' p) («)] = (J {-k\ - \ P X - A}^ (J {ipA - fcA}* (2.34) 

k=0 k=0 

(III) q > p+ 1 : 

zero[t M (u)) = 

poleft^u)] = U {ipA - k\} N (J {-kX + \ P X - qX} N (2.35) 

k=0 k=l 



for the twisted boundary condition and 
(I) q<p-l : 



zero [t^ q \u)\ = set 



i 



pole[t^(n)] = {-qX-\X}{\X} 



2 
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q-l q-l 

\j{-k\-y\-\} 2N \j{y\-k\} 2N (2.36) 

k=0 k=0 



(II) q = p : 
zero[&' p \u)} = {-\p\} 2N 
pole^u)] = {-p\-\\}{\\} 



2 J L2' 

P \J{-k\ - y\ - \} 2N P []{\p\ - k\} 2N (2.37) 

fe=0 k=0 

(III) q > p + 1 : 

zero[t M (u)] = {\ P \-q\} 2N 
pole[*<M)( tt )] = { _ gA -±A}{±A} 

U {ipA - A;A} 27V U {-k\ + ipA - q\} 2N (2.38) 

fc=0 fe=l 

for the open boundary condition, where set g = { — |A} for q = 1 and set g = for q > 1. 

The zeros or poles with order 1 have less contribution than those of order J\f when 
the system size N becomes large. Especially, in the thermodynamic limit iV — > oo only 
these zeros or poles with order N are important. 



3 Functional relations in TV —> oo 

The finite-size corrections for the eigenvalues can be obtained by solving functional 
relations (|2.6| ) and (|2.9|) in the physical strip, 



-A < Re m + \p\ < A . (3.1) 
Denote the finite-size corrections of by T^- tc (u) and write 

r w (t») = rW.(«)T 1 W(Ti). (3.2) 

The bulk and the surface energies determined by the unitary conditions of R and K 
matrices and satisfy 

T^(u)T(±(u + X) = 11 f k (3.3) 

k=0 



Inserting ( |3.2|) into ( |2.27 ) we find that 

rSL(«)ra.(« + A) = i + tW(«). (3.4) 
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So the finite-size corrections for (it) are represented by the hierarchy (it) (t-system 
or i/-system are also called). In the following subsections the analytical treatment of Q3.4Q 
and (|2.9| ) is given. We will see that the finite-size corrections in scaling limit are only 
dependent on the braid asymptotics and bulk behavior of the functional relations. 



3.1 Nonlinear integral equations of real variable 



The Bethe ansatz equations (|2.17 ) render T^ p \u) to be analytic. Since all functions are 



7r-periodic, the analyticity domains for T^'{u) are not unique. It is useful to introduce 
functions of a real variable by restricting the eigenvalue functions to certain lines in the 
complex plane, 

T«(x) := Ti± (%A + P ~l +l x - 3**) , (3-5) 

a*{x) := t {g) (j-x\ + - y>\) and A q {x) : = 1 + a\x) . (3.6) 

The functional relation ( |3.4|) can then be rewritten in terms of the new functions as 

T p (x - lni)T p (x + §7ri) = A p (x) . (3.7) 

For the ground state the functions A^ p \x) and T^ p \x) are analytic, non-zerM in 
— 37r/2 < Im x < 7r/2 and possess constant asymptotics for Re x — > ±oo (ANZC), 
which can be seen from the eigenvalues directly. Taking the logarithmic derivative of the 
above equation and introducing Fourier transforms 

1 POO 

' ■ •• M/ „— ikx 



B p (k) = —j dx [lnT p (x)]' e 

/oo 
dk B p (k) e ikx (3.i 
-oo 



with analogous equations for A g and its Fourier transform A q , then we have 

Transforming back and defining the kernel k(x) 

k(x) := -— — , (3.10) 

v ; 27rcosh(2x) v ; 

we are able to express T q in terms of A 9 , 

lnT p = k * In A p + C p , (3.11) 



1 for those of order TV 
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where C p are integration constants. The convolution / * g of two functions / and g is 
defined by 

poo poo 

(/ *g):= f{x - y)g(y) dy = g(x - y)f(y) dy . (3.12) 



In case of the low-lying excitations states we have to take care of zeros in the analyticity 
strips so that the simple ANZC properties hold. The result (|3.11|) is still correct if we 
change integration path £ so that T p (x) has an ANZC area and Cauchy theorem can be 
applied like the discussion in [Q. The integration constants in ( |3.11| ) can be evaluated 



from the asymptotics of A q and T 9 . In this limit (|3.11| ) becomes 



lnn^fln^ + C*. (3.13) 

It can be seen that the constants are just the multiple of m and do not contribute to the 
jf corrections. 

The A q can be solved from the hierarchy from (|2.9|) , which can be rewritten in terms 
of a q as 

a q {x - \n\)a q {x + ±7ri) = A q - 1 {x)A q+1 {x) . (3.14) 



According to section 2/2 the analyticity strip ( |3.1|) for (u) contains a zero of order M 
at u = — \p\ and a pole of order M at u = — ^p\ + X or u = —^pX — X. All other functions 
are analytic and non-zero in their analyticity strips —\pX — X < u < —\pX + A. We 
introduce finite-size correction terms l q (x) by writing ofl(x) as 



l q (x) , q y£ p 

ta.nh^ (^x)l q (x) , q = p 



"'wh:",,,,.,., vrF (3.15) 



The factor tanh A ^(|x) gives the right zero and poles and all the functions l q (x) therefore 
are ANZC in — 7r < Im x < n. They satisfy the functional equations 

l q {x - \ni)l q {x + ±7ri) = A q -\x)A q+1 {x) . (3.16) 

Again applying Fourier transforms to the logarithmic derivative of the equations and 
then integrating the equations back we obtain the nonlinear integral equations 

lna q = hie 9 + k * ln^ 1 + k * lnA q+1 + D q , (3.17) 

where 

e q {x) := J 1 ' q ^ P (3.18) 

I tanh (^x) , q = p . 

D q are the integral constants. For the same reason we have to take care of the extra 
zeros in the analyticity strips so that the ANZC is held in ( |3.17 ). 
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3.2 Finite-size correction and scaling limit 



The information of finite-size corrections can be extracted from the nonlinear integral 
equations Q3.17 ) and (|3.11|) . The system size N enters the nonlinear equations ( p . 1 7|) 
through (ETTS). The function e p has three asymptotic regimes with transitions in scaling 
regimes when x is of the order of — In iV or In N. We suppose that a q and A q scale 
similarly. So in the following scaling limits, 



e q ±{x) : 

v ' N 



lim e q (±{x + kiN)^ , 
a q ±(x) := lim a q ( ±(x + InM)) , (3.19) 
A q ± {x):= lim A q (±(x + \nAf)) = l + a q ±(x) . 
In this scaling limits, ( 3.17| ) takes the form 

fa q = k q + k * iA^ 1 + k* £A q+1 + D q , (3.20) 
where we use the abbreviations 

h q (x) := lna q {x) , £A q {x) := \nA q {x) , 

*•(»):=( °" "**' (3.21) 

I — 2e , q f 1 p 

and suppress the subscripts ±. The transfer matrix T p (x) in N — ► oo now becomes 

lnT p (x) = {k*\nA p ){x) 

1 r°° ( \nA p (y + \nM) , \nA p (—y — lnJV) 

2lT J- InM 

Mn J-od C ~" ' Nix J~oo~ -'-^v ~" 1 ~\j\f 

Above equation converges and actually can be evaluated explicitly with the help of the 
dilogarithmic function 

L(x) = - [ X dy +|lnxln(l -x). (3.23) 

Jo y 

Multiplying the derivative of (|3T20| ) with £A q , and ( |3T20D itself with (£A q )', taking the 
difference, summing over q, and finally integrating we find 

/oo 
[{h q )'£A q - fa q {£A q )'] dx 
-oo 

/oo 
[{k q )'lA q - {k q - D q ) (£A q )'} dx , (3.24) 
~°° 
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r°° / \nA p {y + \nN) + \nA p {-y - InN) \ ^ | / 1 \ 
J-inAf \cosh(:c — y — \nj\f) cosh(x + y + In J\f) ) \M ) 

/oo p— x /-oo / 1 



where the sum is over all fusion levels q and the contribution of the kernel cancel due to 
the symmetry 



k(—x) = k(x) . 



(3.25) 



Then inserting (|3.21| ) into the right-hand side and integrating the left-hand side of Q3.24| ), 
we are able to obtain 



2 / ^-Wfe) <%=-£!(- 



where the constants D q are given in terms of 



(3.26) 



D q = h q - \HA q ~ x - \PA q+1 



(3.27) 



in asymptotics x —> oo. 



3.3 Asymptotics and bulk behavior 

The nonlinear integral equations (|3.20|) can be easily solved for the limit x — > ±oo and 

A = j- h = 3,A,--- . (3.28) 

For different h it corresponds to different models. The equation ( |3.26| ) shows that these 
asymptotic solutions are enough to obtain the finite-size corrections of the transfer matrix 
T^ p \u). Before discuss the asymptotic solutions it is useful to observe that 

T^-^iioo) =0 for A = n/h , (3.29) 

which can be easily seen from the eigenvalues ( [2.16 ) and the theorem 2.1. 



It is obvious to see that the asymptotics x — > oo corresponds to the braid limit of 
u — > ±ioo. In this limit (12. 9) reduces to 



(^ ) ) 2 = (l+^ 1) )(l + ^ +1) )- (3-30) 
This equation in turn means 

2h q = 0A q ~ l + lA q+1 + D q (3.31) 

in terms of the functions a q . Where the constants D q can be zero or non-zero because the 
different branches can be taken for the dilogarithmic functions in the nonlinear integral 
equations. 
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To solve the let us write tQ as 



sin(3fl) 



sin(0) 

with the parameter # to be determined. The recursion relation ( |3.3U| ) implies 

sm(q9) sin( (q + 2)9 



t 



(?) 



sin 2 # 
sin 2 ((g + l)#) 

*S + 1 = \ 2 ff J (3-32) 

sin v 

for all q = 1, 2 • • -. This solution have to be consistent with the braid limit of T^ q \u). To 
fix the constant parameter # let us consider the "ground" state M = \pN , 



By the relation 



lim T {1) (u)/<j)(u) = 2 cos(y) . (3.33) 



±m = t d) = lim 

Sin(0) °° Im«->±oo 

, 2 /tt\ 

lim -5— ^ 1 = 4 cos 2 - - 1 3.34 



we have 



6 = X = j- . (3.35) 
Moreover the special value of 9 leads the closure condition 

tt~ 2) = . (3.36) 

For the sector j = ^pN — M we have to modify 9 to be 

9 = m\ = ^ (3.37) 
h 

where m = 2j + 1 = 1, 3, • • • < h — 1. For the periodic case m — 1, 2, 3, • • • , h — 1. 

In the limit of x — > — oo can be considered as the bulk behavior in large N. 
According to section |2.2| the analyticity strip for t^ p \u) contains a zero of order N at 
u = —\pX and poles of order N at u = —\pX ± A. All other functions are analytic 
and non-zero in their analyticity strips in —\pX — X < u < \pX + A. For large N the 
leading bulk behavior to the t^ q ' we find that 

constant , q ^ p , 

bulk constant (ta.n(^hu)^j , q — p. ^ 
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The constants are fixed by the functional equations ( |2.9| ) and can be calculated similarly 
to the asymptotics of these t^ 9 '. Like Al model |]4"2" |, it is easy to see that the limit 



lim lim t^' ~ lim exp (— 2e 

x^—oc N—*oo x— >— oo 







(3.39) 



Therefore the functional equations (|279|) are divided into two parts and we find the 
constants for 1 < q < p — 1 



Aq) 

''bulk 



sin(g<j) sinf (g + 2)cr 



sin 2 a 



Aq) , i 

''bulk + 1 



sin 2 ( (g + l)cr 



sin 2 a 



where 



a 



VTL'TX 



p + 2 



m' = 1,2,-- 



Similarly, for p + 1 < q < h — 3 we suppose 



''bulk 



P+ 1 



sin (g — p)r I sin( (g — p + 2)r 



(3.40) 



(3.41) 



sin 2 t 



Aq) , -I 

''bulk + - 1 



sinf (g — p + l)r 



with 



r 



sin 2 t 



m"TC 



h — p 



(3.42) 



(3.43) 



which is consistent with the closure condition (|3.36|) . The eigen-spectra of the transfer 
matrices have only one "quantum" number M, which is related to the braid limit. These 
m' and m" can not be free parameters. In |42| an interpolate method is applied to 
compute the finite size corrections of transfer matrices for ABF models. This follows 
that the exponents m! from the bulk behavior are no longer independent, 



m 



m — m" + 2n + 1 



with the integer n given by 



n 



m — m r 



V 



(3.44) 



(3.45) 



where the brackets \_x\ denotes the greatest integer less than or equal to x. Here for the 
periodic boundary system we need two parameters m, m" and we suppose that m" = 
1,2, ■ ■ ■ , h — p — 1 and m' is given by (3.44). For the largest eigenvalue (or the ground 
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state), the appropriate choices are m = TO' = m" = 1 and m, TO',TO" > 1 give the 
low-lying excited states. The open boundary systems possess the U q [sl{2)] invariance. 



According to the study of XXZ-chain [] 1 we modify ( 3.44 ) to be 

TO — 1 



to' = to + 2n 



n 



P 



(3.46) 



or suppose that to" = 1 and then to' is determined by to. The low-lying excited states 
are given by to, to' > 1. 



The solution t^ lk (u) is given by 

+ A) 



1 < ''bulk )\ L ''bulk J 



(p-ih 



16 cos 2 cr cos 2 r 



Thus we find lastly that 



4mlk( M ) 



1 \AT 

±4 cos cr cos r ( tan(—hu) 



(3.47) 



(3.48) 



3.4 Central charge and conformal weights 



The finite-size corrections are only dependent on the braid and bulk limits of the models. 
In these limits the functional relations are truncated and therefore the sum in ( [3.26D can 
be replaced with the finite sum 



e- y £A p {y) dy 



h-3 , i 



h-3 
9=1 



(3.49) 



where the constants D^'^ can be zero or non-zero because the different branches can 
be taken for the dilogarithmic functions in the nonlinear integral equations. The choice 
of branches have to give the right finite size corrections. Simply taking D^ b ' q ^ = is 
consistent with the asymptotics solutions given in subsection |3.3| . To take nonzero 
we have to single out the right branches of the dilogarithm for the asymptotic solutions 



of the equations, which have been shown for ABF models in [42 



The following useful dilogarithm identity has been established by Kirillov |[45|| . Con- 
sider the functions 

sin(g + 2)(p sin(c/<£>) 



y iq \j,r) :-- 



sin 2 (y9) 



with 



2 + r 



1 < b < n — 1, 1 < q < r 



< j < r 



(3.50) 



(3.51) 
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It is obvious that they are the asymptotic solutions of the functions equations ( |3.20p with 
r = h — 2 or the bulk behavior of the functions equations with r = p and r = h — 2 — p. 
Then the following dilogarithmic function identity holds, 

1 



s(j, r) 



9=1 

L(l 



3r 6j(j + 2) 
+ r 2 + r 



+ 6j) . 



(3.52) 



In terms of the dilogarithm function the finite-size corrections fl3.22| ) are expressed as 



lnT^( 



x 



cosh x ( y, 3 / 1 
Mix V ^ Vl« 



/i-3 
9=1 



(3.53) 



Here we take the case of D q ^ 0. Note that the nonlinear integral equations Q3.20D 
including the closure condition ( |3.36| ) and their solutions presented in subsection [T3] are 



the same as ones of the ABF models studied in |42| . Therefore we can calculate the finite 
size corrections in the same way. Similarly to [EJ, it can be shown that in terms of 
the functions s(j,r) the finite-size corrections ( |3.53| ) for the open boundary systems can 
be expressed as 

it cosh X 



lnT (p) ( 



X 



67V 



s(0,h- 2—p) + s(y,p) - s{m -l,h-2) 

6(1 - m)(p+ 1 - to) + 6u(p -v)\ + fl_ . , , : , . {] 



where v is an unique integer determined by 



v = m — 1 — 



m — 1 
P 



P 



P 



AT 



Inserting ( |3.52|) into ( 3.54j) we have the finite-size correction 



\nT p (x) 



67V ( C 



c — 24A m)J , ) cosh x + o| 



where the center charges c are given by 



3p 



6p 



p + 2 h(h-p) 
and the conformal weights are given by 

2 



(3.55) 



(3.56) 



(3.57) 



(h-(h-p)rnj -p 2 v{j) _ v) 



(3.58) 



Ahp(h - p) 2p(p + 2) 

where p = 1, 2, • • • , h — 2, to = 1,3, ••• < h — 1. Taking into account the geometrical 
factor cosh a; = sin(w/i) and sinhx = icos(hu) we obtain the expression of finite-size 
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correction to the energy given in [ 16] , H7| . For the special case of p = 2 and m — 1 — even 
the same result has been claimed in . 

For the twisted boundary case we have 
7r cosh x 



lnT (p) ( 



x 



(s(m" — 1, h — 2 — p) + s(z/,p) — s(m — 1, h — 2) 
6(m" — m)(p + m" — m) + 6z/(jo — ^ 



— f 

6A/A 



c-24A m „. 



cosh x + o| 



where is an unique integer determined by 



v = m — m" — 



m — m ,j 



P 



P ■ 



(3.59) 



(3.60) 



The conformal weights are given by the following standard expression 



(ht-(h-p)s^ -p 2 v{v _ v) 



4hp(h — p) 



2p(p + 2) 



with 



v = s — t 



s-t 
p 



P 



(3.61) 



(3.62) 



given in |50| , [42]. Comparing with the standard expression we find that t — l,s - 
1, 3, • • • < h— 1 for the fused six- vertex models with the open boundary or t — 1, 2, • • • , h- 
p — l,s — 1, 2, • • • , h — 1 for the fused six- vertex models with the twisted boundary. 



4 Discussion 



We have constructed the functional relations among the fused transfer matrices of the six- 
vertex model with the open boundary conditions. The fusion procedure shows that the 
fusion level of the model can be any positive integers and therefore the fusion hierarchy 
of the six-vertex model with periodic boundary conditions is infinite. The functional 
relations correspond to the su(2) fusion rule without truncation. It means that the 
underlying algebra su(2) has infinite level. This theory has the central charge 3p/ (2 + p). 
By open boundaries, however, the situation can be dramatically changed. The functional 
relations for the six vertex model with the open boundary in braid limit are truncated for 
qh _ e /nA _ g Q su(2) fusion rule with a finite level appears again for the model. 
This shows that the open boundary lower the central charges to be less that 3p/(2 +p). 
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The central charges of underlying conformal field theories for the fused six-vertex 
models with the open boundary condition have been found to be the same as those for the 
fused ABF models. The open boundaries ensure U q [sl(2)} invariance and also play the role 
of the charge at infinite in the Feigin Fuchs construction J57| . The conformal weights of 
the models now take only subset of the Kac formula ( |3.61|) or A Sj „ = Ai jU}S . The quantum 
number s from the braid solutions of the model is specified by the symmetric algebra 
U q [sl(2)]. That means that the spin S enters the calculation through the asymptotics 
of inversion identity hierarchies t q . Suppose that S z = j. Then we have s = 2j + 1 



| 35fl . Specially, we have known that the partition function is a single form in Virasoro 



characters JD], [53], [J5|. 

For the models with the twisted boundary we have found the same central charges 
and the conformal weights. But U q [sl(2)] is no longer the "symmetric group". With the 
similar analysis to [35] we can find that the states of T m {u) with different m are mapped 



on each other by the U q [sl(2)}. We need two integers s and t for the representations of 
Virasoro algebras. As we know the partition function is a sesquilinear form in Virasoro 
characters. 

The six- vertex model with the boundary specified by the reflection matrices K + {u) = 
K(-u - A, f+) and K_(u) = K{u, f_) 

K(u)=( Sm{u + , ° \ (4.63) 
\ sin(-u + £) J V ' 

is the original model studied by Sklyanin ||. This model does not possess U q [sl(2)] 

71 

invariance. In the case of £ + + = — we can solve the finite size corrections similarly 

h 



and thus the similar conformal spectra ( fj.57| ) and ( 3.58| ). One fact to see this is to quickly 



check the braid limit of the transfer matrices, which are the same as those of the model 
with the open boundary condition (|1.15|) . 



The conformal spectra given in this paper is for the fusion hierarchy and therefore 
is more general. For the case of p = 1 the underlying model is the unfused six-vertex 
model. The conformal spectra coincides with that of spin-| XXZ chain with open 
boundary |9|, |35|. The conformal spectra for the general fusion level p gives the spectra 



of spin-ip XXZ chains with the open boundaries. 

The analytic method to calculate the finite size corrections of transfer matrix by solv- 
ing the functional equations has been described for study of the ABF models [j4~2"| |. Here 



we have generalized the method to find the finite size corrections of transfer matrix with 
open boundaries. We like to mention that there is other method available to calculate 
the finite size corrections of transfer matrix with open boundaries, which generalizes the 
method described in pO, IBTI . 
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Appendix A: Fused weights of the six-vertex model 

We give the explicit expression for the fused weights of the six-vertex model in this 
section. Let Y p be the projector on the space of symmetric tensor in C 2p , 

Y = I(pi* + ... + pp-i* + j)...(pi.2 + j) 
pi 

pi,j = ijM(o)/ s in(A) . (A.l) 
The fused weight R p , q (u) G C 2s+1 <g> C 2s+1 is defined by 

R (p,q)( U ) = Y qR(p,q)(u-q\ + \)---R( J> ,2)(u-\)R {pA) (u)Y q (A.2) 

R (pd) (u) = Y p R 1 ' j (u)R 2 ' j (u + X)---R p ' j (u + p\-X)Y p (A.3) 

The derivation of the fused R matrix is straightforward and it is (p + l) 2 by (p + l) 2 
matrix with the following elements 

R ( P ,q)( u )ij = c ( u )J2( F ( n )ifj F ( n i u )ij) *».7>M = s,-s + l.~-,s (A.4) 



25-1 , 25-1 

m=0 v n=l 



C{u) = f[ (sin~ 1 (2sA-mA) \\ sin(w + nA - mA) ) (A.5) 



, , k i sm(s+i+n — m+l)A sin(s— i — m+l)A sin(s— j— m)A / Afi x 

^ = iA " sin(mA) " _ i = \ " sin(mA) " J = \ sin(2 S -m)A ( " } 



, , k i s ^ n sm(u + s\—iX — n\ — m\ + \) s V-^ n sm(u+i\+j\—m\ + \) , , 
(n ' = iA ^M) ii sim>A) (A ' 7) 

where s = and the sum over n is from max(0, j — I) to min(s + j, s — I, s — i). The 
matrix R Pt i(u) can be written as 2 by 2 matrix 

D /x tt • t i \ \ V" f sin(u+p\-m\)e m , m sin(pA-mA)e m4 .i im \ 
R(pi){u) =11 smfw+nA) > (A. 8) 

n= i m =o \ sin(mA)e m _i im sin(u + mA)e m , m y 

where the matrix e^j G End(C 2s+1 ) has only non-zero entry which is 1. 
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Appendix B: Eigenvalue problem of T^ p,l > 

Here we would like to show the Bethe ansatz solutions of T^' 1 '. The transfer matrix 



T M (z) = tr X + (^ 1 /2 )[/(pi) ^ )ir _ ( ^-i/2 )t/ -i i)(z -i ) 
with the monodromy matrix 



U(p,!){z) 



c,2 



L (1,P) 



(1,P) V 



(B.l) 



(B.2) 



is given by the fusion in the "quantum" space. Removing the zeros generated from fusion 
and replacing u by u — \p\ the fused R matrix reads 



' ' ^ • sin(mA)e m _i jm sin(u — ^pX+m\)e mt 



m=0 \ 

which can be rewritten as 

/ zq 8 ' -z^q-S' (q-q-^S 
where the operators 5 Z ,5 ± are generators of C/ g [si(2)] 



sin(pA-mA)e mf i im 



— 2 q° 



q^S^ = q ± S ± S+S--S-S" 



[2S Z 



(B.3) 



(B.4) 



(B.5) 



The matrix ( |B.4|) is just the L-matrix used in |12[]. Therefore the Bethe ansatz solutions 
of T( Pjl )(z) defined in ( P-l| ) have been given exactly in |12| (also see fll3f), which are the 
equations ( |2.16| ) and ( j2.17| ). 



Appendix C: Functional equations from fusion procedure 

In this section we explain how the functional equations of the fused transfer matrices 
come out from fusion procedure. As a simple example, we only show that 



^(1) /Tl(l) 

■*■ 1 1 



if } + f I 

graphically. For this purpose let us represent the R— and K— matrices by 



R (u) 



Kin) 




u. K+(u) 




Therefore we can represent the Yang-Baxter equation (|1.6| ) by 





(c.i; 



(C.2) 



(C.3) 
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the reflection equation ( 1.12| ) by 




























or 
















(C.4) 



and the unitary relation R l,2 (u)R 1,2 (—u) ~ / by 

1 — 

2 — 




(C.5) 



The transfer matrix T(u) with open boundary ( |1.13|) is then represented by 

K+{u+\\) 



K-{u-\\) 



(C.6) 



Let us consider 



K+{u + \\) 



























it 4 


A 














u 











K_(«+±A) 



X_(«-iA) 



(C.7) 



Inserting the identical operator into the position a, b and using the unitary condition 
( p.5| ) and the Yang-Baxter equation ( |U.3| ) (the spectra parameter u is shifted to be 
u — |A in the monodromy matrix ( |1.8|) ), we are able to obtain 



K + (u+%\) 



K+{u+\X) 




u+X 




The fully symmetric and antisymmetric operators are given by the R matrix, 

Y 2 = \(I + P 1 ' 2 ) = D+R X ' 2 (X) 
Y~ = \(I -P 1 ' 2 ) = D-R^(-\) 
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(C.8) 

(C.9) 
(CIO) 



D 



± 



(C.ll) 



where 

/sin- 1 (2A) \ 

±ism- x (A) 

iisin-^A) 

\ sin- 1 (2A)/ 

Then inserting the identical operator into the position a, c of (|C.8|) and using Y 2 + Y~ = I 
and the fusion of the R matrices 



i + A 




(C.12) 



and the K matrices 



K+(u+ 1 A) 



K+{u+\\) 




~ K { ":' 2 \u+\\) 




(C.13) 



K.{u-\\) 



K^ 2) (u-\\) 



The term involved by Y 2 in 
fusion level 2 



< 2) («+iA) 



(C.14) 



gives the transfer matrix with the open boundary of 



(C.15) 

(2) 

which is T . The another term with the antisymmetric projector Y~ , which collapses 
the matrix, is proportional to the identical matrix 



(C.16) 



Therefore we have the equation ( |U.1| ). It can be seen that the proof is correct also for 
the non-diagonal reflection K matrices. So it is very clear that the su(2) fusion rule also 
works for the six- vertex or the eight- vertex models with non-diagonal reflection matrices. 
These will be published elsewhere. 
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